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Abstract 

We apply the classical mapping technique developed recently by Dharma-wardana and Perrot 
for a study of the uniform two-dimensional electron system at arbitrary degeneracy and spin- 
polarization. Pair distribution functions, structure factors, the Helmhotz free energy, and the 
compressibility are calculated for a wide range of parameters. It is shown that at low tempera- 
tures T/Tp < 0.1, Tp being the Fermi temperature, our results almost reduce to those of zero- 
temperature analyses. In the region T/Tp > 1, the finite temperature effects become considerable 
at high densities for all spin-polarizations. We find that, in our approximation without bridge func- 
tions, the finite temperature electron system in two dimensions remains to be paramagnetic fluid 
until the Wigner crystallization density. Our results are compared with those of three-dimensional 
system and indicated are the similarities in temperature, spin-polarization, and density dependen- 
cies of many physical properties. 

PACS numbers: 73.21.Fg, 71.45. Gm, 71.10.-w 
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I. INTRODUCTION 



Two-dimensional electron systems have been investigated by a number of investigators. 
Their importance seems to have recently increased due to possible applications to spintronics 
for electronic devices. The properties of two-dimensional electron systems are determined 
by the Coulomb interaction between electrons and depend crucially on the density, temper- 
ature and spin-polarization. In principle, exact results can be obtained by the numerical 
experiments and are available to some extent at zero temperature. At finite temperatures, 
however, reliable numerical experiments are still awaited and the importance of theoretical 
approach is not reduced. 

Recently, Dharma-wardana and Perrot have developed a method called classical-map 
hypernetted-chain method (CHNC)^*^^. They determine the quantum temperature Tg so 
that the classical system at the temperature Tg has the same correlation energy as the 
quantum system at T = and assume that the properties of the system at the finite 
temperature T is given by those of classical system at the temperature 



One can include the higher-order cluster interactions in classical fluids applying the modified 
HNC method^ and taking the bridge function into account. 

Various properties of two-dimensional electronic systems at zero-temperature have been 
examined by Bulutay and Tanatar^ based on the CHNC method without bridge corrections. 
We here extend the analysis to finite temperature two-dimensional electronic systems with 
arbitrary spin-polarization. Confirming that our results reduce to those of Ref. 6 at low tem- 
peratures such that T/Tp < 0.1, we analyze the finite temperature effects for T/Tp > 0.5, 
Tp being the Fermi temperature. We observe that the effects are substantial at high densities 
irrespective of spin-polarization. It is shown that, within our analyses without bridge func- 
tions, the two-dimensional electron system remains paramagnetic until the Wigner lattice 
formation. In Ref. 4, the existence of a ferromagnetic phase before the formation of Wigner 
lattice is predicted and, in such a domain, one needs careful treatment including the bridge 
function. The purpose of this paper is to analyze fundamental quantities in a wide domain 
of density, temperature, and spin-polarization where one may neglect the bridge function as 
in Ref. 6. In what follows, we adopt the atomic units and take = 1- 




(1) 
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II. CHNC METHOD APPLIED TO TWO-DIMENSIONAL ELECTRON SYSTEM 



We consider a two-dimensional fluid of electrons at the temperature T containing two 
spin species of surface densities rti and n2 {ni > with the total density n = rii + n2 
and spin-polarization C, = (ni — n2)/n. We determine the Pauli potential for parallel-spin 
electrons Pui^r) by the relation 

glir) = exp[-f3Pu{r) + hl{r) - 4(r)]. (2) 

Here P is determined by the temperature T^f by = l/T^f, fi'°(r) = (r) + 1 is the ideal 
gas pair distribution function for the component i, 

hUr) = E nihHh) exp[^(ki - k^) ■ r] = -[/.(r)]^ (3) 

^« ki,k2 

and c° (r) is the direct correlation function corresponding to (r). In (3), n{k) is the Fermi 
occupation number at temperature T. At T = 0, /i(r) = 2Ji{kpr)/kpr, where Ji{x) is the 
Bessel function and kp is the Fermi wave number of species i. 

The pair distribution functions (PDF's) and the direct correlation functions are related 

via 

9ij (r) = exp [-P (pij (r ) + hij (r) - (r ) + (r)] . (4) 
Here 4>ij{r) is the pair potential between the species i and j defined by 

= Pii{r)6ij + Vconir), (5) 

and hij{r) = gij{r) — 1, Cjj(r), and Bij{r) are the pair correlation, direct correlation and the 
bridge function, respectively. In the pair potential, Vcou{r) includes the diffraction effect as 

Vcouir) = i[l-e-'-^-], (6) 

where kt^ = (27rm*Tc/)^/^, m* = 1/2 is the reduced mass of the scattering electron paii^. 
We also have the Ornstein-Zernike relations between the pair correlation functions and the 
direct correlation functions: 

hijir) = c,j{r) + Y.^s [ dr' /ii,(|r - r'|)c,j(r'). (7) 

s •' 

The equations (jH) and ((Tj) are exact but are not closed because of the existence of the 
unknown bridge functions -Bjj(r). In this paper, however, we neglect the bridge functions as 
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in Ref. 6 and solve the set of equations (0} and (|7j). We have analyzed^ the case of = 5 
with the bridge function employing the solution of the Percus-Yevick equation for hard disks 
with somewhat different formula for Tq. The result indicates that the effect of the neglect 
of bridge function increases with the temperature and the free energy due to exchange- 
correlation is overestimated by at most 5% at T jTp = 5. The results for fundamental 
quantities obtained without bridge functions will therefore be still useful in a wide range 
of density, temperature, and spin-polarization where the effect of the neglect of the bridge 
function is not so subtle. 

In the CHNC method, the temperature of the classical fluid Tcf is determined by ((H). 
For the quantum temperature Tq, we adopt the one given by Ref. 6, 

^ _ 1 + avs riiTpi + n2TF2 
h + CTg n 

where r, = {irn)-^/'^ with a = 1.470342, b = 6.099404, and c = 0.476465. 

The exchange-correlation part of the Helmholtz free energy is calculated by the integration 
over the coupling A as 



xc 

— = nn 
n 



£ dX J dr[g{X;r)-l] (9) 



from the spin-averaged PDF given by 

9{r) = ^[(1 + C)'^7ii(r) + 2(1-C'ki2(r) + (1-Cfe(r)]. (10) 
The total Helmholtz free energy Ftot is obtained as^ 

Ftot = Fo + F^ci (11) 
where Fq = Fq + Fq is the free energy of the ideal electron gas given by 

F^ = n,^^, - El (12) 

/z^ = Tln[e^>/^-l], (13) 

£;i. = (l±C)vrn, (14) 

and 

° 71 Jo exp[x- ni/T] + r ^ ^ 

We calculate the exchange part of the Helmholtz free enegy F^ = F^ + F^ by the ideal 

gas values as 

^ = ^^^1 + 0'/ rfr[^?,(r)-l], (16) 



n 4 J 



(17) 



and define the correlation part by F^. = F^c ~ F^- At T = we have 



n n o-KTs 



(18) 




(19) 



III. RESULTS AND DISCUSSIONS 

We have solved the coupled equations (jH) and ((7j) using the Hankel transform as in Refs. 
6 and 9 and neglecting the bridge functions for a wide range of density, spin-polarization 
and temperature, 0<rs<30,0<C<l, and 0<t = T/Tp < 5. Here we summarize the 
results. 

A. Pair distribution function and structure factor 

The spin-averaged PDF g{r) for = 1 and T/Tp = 0, 1 in two cases C = 0, 1 is shown 
in Fig.l. Our results for T/Tp = are in good agreement with those given in Ref. 6. We 
observe that the finite temperature effect considerably reduces the exchange-correlation hole 
around electrons in both cases of C = and ( = 1 and we may estimate the behavior of the 
pair distribution function in the case of arbitrary spin polarization from these results. 

In Fig. 2, we show the spin-averaged structure factor defined by 



for = 1, T/Tp = 0, 1, and ( = 0,1. At T/Tp = 0, the structure factor is consistent 
with the one given in Ref. 6. We observe the finite temperature effect for the cases of 
spin-polarization = 0, 1. 

B. Total free energy 

We have calculated the total Helmholtz free energy for finite temperatures < T/Tp < 5 
and spin-polarizations < C ^ 1- The results are shown as a function of the density 
parameter and spin-polarization ( in Figs. 3 and 4, respectively. 




(20) 
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At low temperatures T /Tp < 0.1, our results reduce to those of Ref. 6. When T/Tp > 1, 
the finite temperature effects become considerable at high densities for all cases of spin- 
polarization. Our results indicate that the two-dimensional electron system remains to 
be in the paramagnetic fluid phase until the Wigner crystallization densities even at high 
temperatures. 

The dashed-dotted line in Fig. 4 represents the total free energy for T/Tp = 1 and 
Vs = 1. It is seen from the figure that the spin-polarization effect is more pronounced at 
higher temperatures and densities. 

C. Exchange-correlation free energy 

To compare the electron correlation effects in two and three dimensions, we plot the 
exchange-correlation and correlation free energies as a function of temperature in Figs. 5 
and 6. We find that our results of the temperature dependence of the free energy in two 
dimensions are similar to those in three dimensions^. 

To check the validity of the CHNC scheme used in this paper, we compare our exchange- 
correlation energy with Monte Carlo (MC) and Singwi-Tosi-Land-Sjolander scheme (STLS) 
results given in Refs. 10, 11, and 12. The exchange-correlation energies computed by 
different schemes are listed in Table I. We observe that the difference between our results 
and STLS results is considerable (about 10%) only at high temperatures (T > Tp). This 
may be because we have used the expression for quantum temperature Tg fitted to MC data 
at T = 0. To obtain better agreements with STLS results for T > Tp may need to 
include the bridge term and use the results of finite temperature STLS scheme given in Ref. 
12 for fitting procedure. 

As for the exchange-correlation free energy at high densities, the agreement with MC 
and other values can be improved by interpolating Tq more accurately. For example, the set 
a = 24.936526, b = 87.221405, and c = 8.455547 gives better fitting for Tg as shown in Fig. 
7 and better agreement as listed on the last line in Table I. 
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TABLE I: Exchange-correlation energy per electron (in a.u.) of unpolarized 2DEG, —Fxc- For 
comparison, the MC results of Tanatar and Ceperley^'^ and the STLS results of Jonson^^ and 
Schweng and Bohnti^ are given (denoted by a, b, and c, respectively). The results for fully polarized 
2DEG are shown in brackets. The last line is the values based on improved interpolation for Tq. 
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0.693 
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0.685 


0.585 


0.633 


0.359 


0.399 



D. Compressibility 

The compressibility is obtained by density derivatives of the free energy. In order to 
check the accuracy of our finite temperature scheme, we have calculated the compressibility 
as a function of density parameter or spin-polarizations C, for different temperatures. We 
note that at low temperatures our results reduce to those of Ref. 6. Our results shown in 
Fig. 8 indicate that the temperature effect on the compressibility is remarkable only at low 
densities and high temperatures such that T/Tp > 0.3. 

The spin-polarization dependence of the compressibility is displayed in Fig. 9. We observe 
that the effect of spin-polarization increases with the temperature and electron density. 

IV. CONCLUSIONS 

Applying the recently proposed classical-map hypernetted-chain (CHNC) method to the 
two-dimensional electron system, we have calculated the pair distribution function, structure 
factor, Helmhotz free energy, and compressibility at finite temperatures for a wide range of 
spin-polarization and density. It is shown that at low temperatures T/Tf < 0.1, our results 
are almost identical to those of Ref. 6 at T = 0. Our results indicate that correlation char- 
acteristics in 2D electron liquids depend remarkably on the temperature, spin-polarization, 
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and density, similarly to the three-dimensional case. 

When T > Tp, the finite temperature effects become considerable at high densities irre- 
spective of spin-polarization as concretely shown in various quantities in this paper. We find 
that the finite temperature 2D electron system remains to be in the paramagnetic fiuid phase 
until the Wigner crystallization density is attained. Though the validity of the conclusion 
on the phase around ~ 30 might depend on that of CHNC applied to two-dimensional 
electron system and the neglect of bridge functions, results obtained here for 1 < < 30 will 
be useful in investigation of the system whose finite temperature properties are important 
in applications but not exactly known. 
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FIG. 2: Structure factor S{k) for = 1, t = T/Tp = 0, 1, and C = 0, 1. 
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FIG. 3: Total Helmholtz free energy as a function of for different temperatures and spin- 
polarization. 
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FIG. 4: Total Helmholtz free energy as a function of spin-polarization ^ for different temperatures 
and densities. 
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FIG. 5: Exchange-correlation free energy per electron in units of /j^x = — 2A;F/vr as a function of 
temperature for different spin-polarization and densities. 




T/T. 



FIG. 6: Correlation free energy per electron in units of Hx = — 2/ej?/7r as a function of temperature 
for different spin-polarization and densities. 
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FIG. 7: Improved fitting for Tq in comparison witli tfie one given in Ref. 6. 




FIG. 8: Inverse compressibility of unpolarized {Q = 0) and fully polarized (C = 1) phases normalized 
by 2D free Fermion value Kq for t = T/Tp = 0, 0.1, 0.3, 0.5. 
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FIG. 9: Compressibility normalized by paramagnetic value = 0) as a function of spin-polarization 
for t = T/Tf = 0.5, 1, and = 1, 5, 10. 
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